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Abstract 

The present study explores the potential of directly linking phylogenetic identities obtained by cloning and 
sequencing of ITS sequences to dominant ribotypes of molecular community fingerprints to give further insight into 
dominant members of the communities in three Irish grassland soils. The ten most abundant bacterial ribotypes of 
untreated bare soils of three grassland microcosms were chosen to represent the "baseline community" of the 
respective soil. Identities on phylum and order level were assigned to these ribotypes on a weighted basis, by 
matching sequence homologies of cloned ITS sequences with ribotypes of the same fragment lengths +5 bp. 
Results showed that ribotypes were represented by the phyla Acidobacteria, Actinobacteria, Proteobacteria, and 
Firmicutes and the distribution of the ribotype and phylotype communities was shown to be highly site-specific. 
Furthermore the response of dominant bacterial phylotypes to plant species composition, fertilisation and Lolium 
perenne ingression was investigated within a larger microcosm study (Microb Ecol 63:509-521). 
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Main text 

Microbial communities are likely to be influenced by a 
wide range of factors, including environmental and an- 
thropogenic factors. In order to examine those influences 
on the community structure, molecular fingerprinting 
techniques, such as terminal restriction fragment length 
polymorphism (TRFLP) or automated ribosomal intergenic 
spacer analysis (ARISA), are commonly employed (Thies 
2007). While these automated fingerprinting techniques 
provide insight into changes of community structure at 
high resolutions, they give little or no information on the 
community members involved in the change. To tackle the 
question, whether or not changes in ribotype communities 
actually reflect significant changes on a phylogenetic or 
functional level, some studies employ amplification of 
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functional or group-specific genes prior to fingerprint- 
ing (Briones et al. 2002; O'Callaghan et al. 2008), while 
other studies complement fingerprinting with cloning 
and sequencing analysis to identify species present in 
samples (Dickie & Fitzjohn 2007; Shi and Bending 
2007; Shi et al. 2011). This study explored the latter ap- 
proach and its potential of directly linking cloning and 
sequencing data to ARISA-derived ribotypes to provide 
further insight into community changes induced by agri- 
cultural practices in three different Irish grassland soils 
(Liliensiek et al. 2012; for soil properties see Additional 
file 1: Table SI), including a near-neutral grassland 
(hereafter called Burren soil), a mesotrophic grassland 
(Ardgillan soil), and an acidic upland grassland (Wicklow 
soil). The phylogenetic identities of ARISA-derived dom- 
inant members were projected onto the fingerprinting 
data of a larger microcosm study (Liliensiek et al. 2012; 
Additional file 2: Table S2) to investigate the response of 
dominant bacterial phylotypes to plant species compos- 
ition, fertilisation and Lolium perenne ingression. 
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In order to assign phylogenetic identities to the 10 
most abundant bacterial ARISA ribotypes (hereafter 
referred to as baseline ribotypes) of each untreated 
bare soil type (indicated by shaded cells in Additional 
file 2: Table S2), PCR products (primers ITSF/ITSReub, 
Cardinale et al. 2004) of DNA extracted from micro- 
cosms of each of the untreated bare soil types (Liliensiek 
et al. 2012) were randomly cloned (pGEM-T easy, 
Promega), sequenced (primers SP6/T7, Macrogen Inc.) 
and matched to homologies in the NCBI database 
using BLAST (Altschul et al. 1997; see Additional file 3: 
Methods). Clones were then assigned to the baseline 
ribotypes based on corresponding sequence lengths in 
order to confer phylogenetic identities to ribotypes. Se- 
quences targeting the ITS region have previously been 
applied to phylogenetic analyses (Ranjard et al. 2000), 
but contain only short stretches of the 16S and 23S genes. 
Furthermore, sequences in the database are based on 
culturable organisms with potentially low sequence hom- 
ologies to environmental 16S genes (van Wintzingerode 
et al. 1997). Thus not all ribotypes could be exactly 
matched with a phylogenetically identified clone. Tolerat- 
ing a mismatch in length of ±5 bp between ribotype 
number and cloned sequence length, however, about 67% 
of ribotypes present in the baseline-communities could 
be assigned down to phylum level (Additional fUe 4: 
Table S3). Given a certain degree of inaccuracy during frag- 
ment analysis as well as sequencing, such a variation was 
deemed possible to occur. In order to limit the resulting 
bias, matches between ribotypes and their respective clones 
were assigned on a weighted basis. Weighted abundances 
of ribotypes associated to the same phylum were then aver- 
aged to give an estimation of the phylogenetic distribution 
(see Additional file 3: Methods; Table 1). 

Ribotypes present in the baseline-communities were 
represented by the phyla Acidobacteria, Actinobacteria, 
a-Proteobacteria, y-Proteobacteria and Firmicutes, all of 
which have frequently been found in grassland soil (van 
Elsas et al. 2002). The distribution of these phyla among 
the top 10 ribotypes was however, significantly dif- 
fered between soil types (Table 1), as shown by ANOVA 
(Additional file 3: Methods). Actinobacteria and Firmicutes 
were the most common and abundant phyla considering 
all soils combined. These phyla generally include the gen- 
era Streptomyces, Clostridium and Bacillus, which have 
been reported to exhibit the highest copy numbers (10-12 
copies) of the RNA-operon (Tourova 2003). Thus, unless 
ribotypes can be identified down to genera-level, it remains 
unclear if individuals represented by abundant ribotypes 
are indeed highly abundant, or just appear as such due to a 
higher RNA-operon copy number. 

While Actinobacteria and Firmicutes were represented 
by baseline-ribotypes from all three soils, Proteobacteria 
were solely present in the baseline-community of the 



Burren soil which is in agreement with their prefer- 
ence of a more alkaline pH (Kumar and Nicholas 
1983). Acidobacteria, which are moderately acidophilic 
(Hugenholtz et al. 1998), on the other hand were found in 
the baseline-community of the Wicklow soil only. 

Despite their presence in the baseline-communities of 
all three soils, the distribution of Actinobacteria and 
Firmicutes was significantiy affected by soil type, whereby 
all soil properties (sand, silt, clay, pH and nitrogen and 
carbon content) contributed significantly to this effect for 
Firmicutes, but not for Actinobacteria, where pH and its 
interaction with clay were the most significantly contribut- 
ing factors. The reputed preference of Actinobacteria for 
an alkaline pH (Atlas and Bartha 1997) was reflected by 
a significant decrease of their relative abundance with 
pH, from the Burren to the Wicklow soil (pH decrease 
from 6.45 to 4.45 in the order of Burren soil > Ardgillan 
soil > Wicklow soil). These observations agree with pre- 
vious findings that, regardless of soil location, molecular 
studies seem to detect the same major phyla while the 
community structure at phylum level is highly site spe- 
cific (Hackl et al. 2004). 

Applying the phylotypic identities of baseline-ribotypes 
to their occurrence in fingerprints of the complete 
microcosm study (Liliensiek et al. 2012), treatment (i.e. 
plant species composition, fertilisation and presence of 
Lolium perenne; Additional file 2: Table S2) effects on the 
distribution of the baseline-communities were tested. 
Ribotypes defined as baseline communities in untreated 
bare soils generally were among the top 10, sometimes 
20 most abundant ribotypes detected under any treat- 
ment condition within the same soil, their relative abun- 
dance distribution however was affected by treatment. 
Plant species was the main factor influencing most base- 
line ribo- and phyloypes, followed by fertilisation and the 
plant-fertilisation interaction (Table 2). The most pro- 
nounced plant species effect was observed in the Burren 
baseline-community with an influence on all ribotypes 
and phylotypes (F = 53.9 to 154.6 and F = 62.2 to 522. re- 
spectively), except one ribotype (F = 2.5). In the Ardgillan 
community 6 out of 10 ribotypes (F = 5.2 to 9.7) and only 
Actinobacteria (F = 4.1) displayed a significant response 
to plant species, while 5 out of 10 ribotypes (F = 3.9 
to 32.7) and Acidobacteria and Firmicutes (F = 3.7 and 
15.6, respectively) were influenced in the Wicklow 
soil. The Lolium effect was generally very weak, with only 
one ribotype influenced in the Burren and Ardgillan site 
(F = 4.7 and 3.0 respectively) while phylotypes seemed 
not to be affected. The Wicklow community showed no 
response to Lolium ingression, neither on the ribotype, 
nor phylotype level. A fertilisation effect on the phylo- 
type level was only observed in the Burren and Wicklow 
soils, where it influenced Firmicutes (F = 12.3) only, 
and Firmicutes (F = 26.1) and Actinobacteria (F = 10.4), 
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Table 1 Means and variations of relative abundances (rfu) of top 10 ribotypes and associated phyla as well as ANOVA 
F-values for their distribution across bare untreated soils 



Ribotype NCBI matches 


Top 10 of site 




Relative abundances (rfu) 




ANOVA 








Burren 


Ardgillan 


Wicklow 




225 


ni 


W,A,B 


1480±130° 


2704±259 " 


2265±432 ° 


4.26 ™ 


241 


ni 


W 


0±0 " 


538±83 


1796±212 ° 


49.08" 


282 


An, aP, yP 


B 


2982±272 ° 


488±116 


0±0 


87.99" 


301 


An, F 


A 


872±124° 


1444±165 ° 


803±341 ° 


2.34 


302 


An, F 


A 


0±0 


1499±82° 


0±0 


333.49" 


324 


ni 


B 


1323±102 ° 


422±38 ^ 


0±0^ 


1 1 5.58" 


352 


ni 


W 


0±0 


0±0 


1 643±206 ° 


63.84" 


390 


An, F 


W 


0±0 


0±0 


4387±308 ° 


202.24" 


391 


An, F 


W 


0±0 


0±0 


1 634±204 ° 


64.24" 


398 


An, F 


W 


0±0 " 


705±16 


1538±328 ° 


1 6.49"* 


400 


An 


A 


647±48 


3704±106° 


1158±237 '' 


1 1 5.20"** 


409 


An 


W 


0±0 " 


853±108 ^ 


1667±29° 


168.33*** 


417 


An 


B 


1 938±63 " 


0±0 " 


214±65 


411.27*** 


502 


yp 


B 


1427±35 " 


0±0 


0±0 ^ 


1702.14*** 


612 


ni 


B 


4103±85 ° 


0±0 


0±0 


2330.72*** 


831 


F 


A 


885±96 


2327±132° 


0±0^ 


1 54.64*** 


835 


F 


A 


1 1 90±83 


2002±175 ° 


0±0^ 


80.81*** 


843 


F 


A 


0±0 ^ 


1183±76° 


0±0 


244.56*** 


846 


ni 


A 


0±0 


3296±153 ° 


0±0 ^ 


463.38*** 


850 


F 


B 


4376±476 ° 


766±122 


0±0 ^ 


67.94*** 


853 


F 


B 


3241 ±21 2° 


0±0 


0±0 


234.04*** 


890 


ni 


A,B 


2742±125 " 


2021 ±51 2° 


0±0 


21.82** 


895 


F 


B 


1411+121 " 


0±0 


0±0 


135.38*** 


904 


F 


A 


953±13 


1167±95 " 


0±0^ 


126.18*** 


917 


ni 


W 


509±121 


0±0 " 


2216±118° 


141.19*** 


920 


ni 


W 


0±0 


0±0 


9694±612 ° 


250.79*** 


953 


Ad 


W 


0±0 


0±0 


1555±258 ° 


36.33*** 


Phylotype 




Acidobacteria 


W,A,B 


0±0 


0±0 


1 3998±2322 ° 


36.33*** 




Actinobacteria 


W 


1 7202±678 ° 


1 1 368±240 


9051 ±500 " 


69.05*** 




a-Proteobacteria 


B 


16234±1365 " 


0±0 


0±0 


14136*** 




Y-Proteo bacteria 


B 


17702±1239 ° 


0±0 


0±0 


204.12*** 




Firmicutes 


W,A,B 


13325±384 


8887±191 " 


17157±890 ° 


52.66*** 



Means and variations of the relative abundances in relative fluorescent units (rfu) of individual ribotypes and associated phyla (see Additional file 3: Methods) are 
presented for each site. Column three indicates the membership in the respective baseline community (Burren (B), Ardgillan (A), Wicklov\/ (W)}. Superscript letters 
indicate the group differences between relative abundances across sites. The significance of distribution differences between sites are presented in the last 
column by F-values and their significance level (p = 0.05 > * > 0.01 > > 0.001 > ***). 

Abbreviations: ni - not identified. An - Actinobacteria, Ad - Acidobacteria, F - Firmicutes, aP - a-Proteobacteria, yP - y-Proteobacteria. 



respectively. Similarly 5 and 6 ribotypes were influenced in 
these soils, respectively (F = 4.2 to 11.8 and 4.8 to 30.0, re- 
spectively), while only 1 ribotype was affected by fertilisa- 
tion in the Ardgillan soil (F = 15.8). 

In summary, the constitution of the baseline community 
on ribotype and phylogenetic level showed some unique 



but mostly common members of each site and the distri- 
bution of relative abundances on both levels was signifi- 
candy affected by site and treatment. While Proteobacteria 
were found solely in the baseline-community of the near 
neutral pH Burren soil, Acidobacteria were only found in 
the acidic Wicklow soil. The baseline-ribotypes in the 
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Table 2 ANOVA table of F-values for ribotype and phylotype responses to plant species, fertilisation and Lolium 
perenne ingression as well as their interactive effects 





Ribotype 


Phylotype 


Plant 


Fert 


Plant X Fert 


Lp 


Plant X Lp 


Fert X Lp 


Plant X Fert x Lp 


Burren 


225 


ni 


2.5 


1.8 


1.3 


0.0 


4.2'* 


0.9 


1.7 




282 


An, aP, yP 


53.9™ 


1.1 


0.7 


0.0 


1.0 


2.1 


2.3 




324 


ni 


181.3™ 


0.3 


3.9" 


04 


24 


0.6 


0.1 




417 


An 


88.7™ 


0.9 


1.5 


0.7 


2.1 


4.5* 


4,3** 




502 


yp 


1 07.2™ 


5.5* 


1.0 


0.3 


1.5 


0.2 


0.5 




612 


ni 


1 1 8.4™ 


1 1 .7" 


3.2* 


4.7* 


1.1 


0.2 


0.4 




850 


F 


140.2™ 


1 1 .8" 


5.5" 


1.0 


1.0 


4.2* 


1.2 




853 


F 


1 29.6™ 


1 1 .0" 


6.8*" 


0.1 


3.2* 


0.0 


1.1 




890 


ni 


1 54.6™ 


0.7 


0.9 


1.0 


1.7 


0.2 


0.5 




895 


F 


1 34.3™ 


4.2* 


4.3" 


1.8 


7.r* 


0.9 


3,0* 




Actinobacteria 


1 33.0™ 


0.1 


14 


0.1 


2.3 


0.0 


0.5 




Firmicutes 


522.0™ 


1 2.3" 


9.3*" 


0.1 


3.0* 


04 


05 




a-Proteobacteria 


62.2™ 


1.0 


0.8 


0.0 


1.1 


1.9 


2.2 




y-Proteobacteria 


92.6™ 


2.5 


0.6 


0.0 


14 


1.9 


2.1 


Ardgillan 


225 


ni 


6.3™ 


2.1 


7.7*" 


0.0 


0.8 


0.3 


2.2 




301 


An, F 


8.4™ 


0.3 


6.2*" 


0.6 


5.6'* 


5.7* 


6.9*** 




302 


An, F 


1.7 


0.5 


0.3 


0.5 


0.9 


0.8 


0.6 




400 


An 


5.2*' 


0.0 


2.1 


0.6 


0.1 


04 


1.0 




831 


F 


8.7™ 


04 


0.7 


3.1 


0.7 


2.9 


0.6 




835 


F 


8.6™ 


0.1 


0.6 


1.3 


0.7 


34 


0.2 




843 


F 


2.1 


0.2 


1.0 


1.5 


1.2 


0.3 


0.7 




846 


ni 


9.7™ 


15.8*" 


3.0* 


7.6" 


1.7 


0.8 


0.3 




890 


ni 


2.5 


2.2 


04 


0.1 


1.1 


0.3 


0.3 




904 


F 


14 


04 


0.9 


0.2 


2.0 


0.1 


05 




Actinobacteria 


4.1" 


0.9 


6.5*" 


0.2 


1.3 


1.1 


2.8* 




Firmicutes 


1.3 


0.6 


4.2" 


0.0 


1.0 


2.5 


0.6 


Wicklow 


225 


ni 


1 3.8™ 


6.9* 


9.8*" 


0.3 


1.1 


3.1 


3.9" 




241 


ni 


3.9" 


7.1* 


1.6 


1.2 


1.8 


0.0 


1.0 




352 


ni 


32.7™ 


30.0*" 


9.9*" 


0.1 


1.6 


0.1 


2.4 




390 


An, F 


1 0.8™ 


12.9™ 


1.9 


0.9 


1.0 


9.r* 


3,8* 




391 


An, F 


24 


4.8* 


0.5 


0.6 


1.0 


0.0 


1.6 




398 


An, F 


0.5 


4.1 


0.5 


0.1 


1.0 


2.3 


1.7 




409 


An 


2.3 


0.3 


04 


0.1 


1.0 


0.0 


1.3 




917 


ni 


1 2.0™ 


9.2" 


7.1*" 


0.3 


2.3 


0.3 


6.5"* 




920 


ni 


1.9 


2.4 


2.0 


1.0 


2.3 


2.6 


0.4 




953 


Ad 


250227.0 


52097.1 


190238.3 


1 7340.0 


81344 


49192.1 


1083204 




Acidobacteria 


3.7* 


0.8 


2.9* 


0.3 


0.1 


0.7 


1.6 




Actinobacteria 


24 


10.4" 


1.2 


0.8 


1.0 


4.0 


1.5 




Firmicutes 


1 5.6*" 


26.1*" 


2.0 


2.0 


1.1 


8.7" 


2.3 



Analysis of treatment effects on the baseline-connmunities of respective sites covered the full model of plant species (Plant), fertilisation (Fert) and Lolium 
ingression (Lp), including all their interactions (x). Ribo- and phylotypes belonging to the baseline-community of each site are listed in the first column, followed 
by F-values for each model term and their significance (p = 0.05 > * > 0.01 > ** > 0.001 > ***). 

Abbreviations: ni - not identified. An - Actinobacteria, Ad - Acidobacteria, F - Firmicutes, aP - a-Proteobacteria, \? - y-Proteobacteria. 
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Wicklow soil were less responsive to plant species and 
more susceptible to fertilisation compared to the other 
sites, while plant species, including Lolium, seemed to have 
the strongest effect in the Burren soil. In general interactive 
effects were less strong and common, suggesting a greater 
response to individual treatments. Critical evaluation 
brings up several issues associated to such an approach. 
Obviously, the specificity of identification is low due to 
limited homologies for environmental ITS genes in the 
database as well as the target sequence of ARISA clones. 
More interestingly, communities were dominated by 
Actinobacteria and Firmicutes which include genera with 
the highest RNA-operon copy numbers, raising the ques- 
tion, if studies considering only abundant ribotypes might 
overlook actually abundant species of the community. 
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